%
% interp_matrix_v.m
%
%   Form the interpolation matrix for the vertical velocity (v)
%   which is stored at the horizontal edges
%
%   input -- Nx number of cells in x-direction, fine grid
%            Ny number of cells in the y-direction, fine grid
%   output
%
function Iv = interp_matrix_v(Nx,Ny);

% form two one-dimensional interpolation matricies
%
Ipx = cc_interp1(Nx);
Ipy = ec_interp1(Ny);

% assemble the matrix for the two-dimensional problem
%   note the y-x ordering in kron, this correpsonds to 
%   the case when the first index is for the x-coordinate
%
Iv = kron(Ipx,Ipy);

